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Abstract 

Statistics of Poincare recurrence for a class of circle maps, including 
sub-critical, critical, and super-critical cases, are studied. It is shown 
how the topological differences in the various types of the dynamics 
are manifested in the statistics of the return times. 
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1 



1 Introduction 



Statistics of the Poincare recurrences, i.e. return times statistics, has recently 
gained renewed importance in the theory of dynamical systems, primarily due 
to the fact that it could be used as an indicator of the statistical properties 
of the system's global dynamics on large parts of its phase space (see for 
example PP,|2])- For instance, the first return time can be used to calculate 
the metric entropy of a system with an ergodic invariant measure [3], and 
seems to be related to other generalized dimensions used to describe the 
fractal properties of the dynamics, at least for certain types of dynamical 
systems jjj. In this paper we study the Poincare recurrences for maps of 
the circle, that can display quasi-periodic, bi-stable or chaotic dynamics, 
depending on the values of their two parameters. 

Given a discrete dynamical system on a phase space M with a transfor- 
mation T : M — > M and a reference measure \x on M, the first return time, 
in a measurable set A C M, of a point x G A, is defined by 

t a (x) = mf{T k (x) e A}. (1) 

k 

The first return time ta for the set A, and the average return time < ta > 
for the set A, are the following 

t a = inf [r A (x), <r A >= T A {x)d^ A {x), (2) 

x£A J A 

where fJ, A (x) is the conditional measure ^a(B) = f/,(B) / f/,(A) for any B C A. 

With < ta > and ta(x) we define the probability Fj\,(t) that the normal- 
ized first return time in the set A is larger than t for the points in the set 
A 

F A (t) = fi A (A >t ) A >t = {x e A t : t a (x)/(t a ) > t} . (3) 

A limit probability measure F x (t) may be associated to any point x by con- 
sidering a sequence of neighborhoods A k of x with n(A k ) — > as k — > oo. 
For ergodic systems the limit is the the same for almost every point and in 
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this case the average return time in a domain A is given by 

according to the well known Kac's lemma. For some classes of hyperbolic dy- 
namical systems |S] it has been proved that the return times spectrum follows 
the exponential-one decay law F x (t) = e~* at almost every x. If x belongs to 
the dense set of the unstable periodic points then F x (t) = a exp(-at) where 
a and a depend on the period. It is known that the properties of the distri- 
bution of the return times can give criteria for the existence of an equilibrium 
and the rates of mixing 7\. A polynomial decay law for F x (t) was found for 
integrable area-preserving maps [EHH1- O n the basis of numerical computa- 
tions, the polynomial decay of return times spectra F x (t) was also claimed 
to be a generic property of Hamiltonian systems with mixed phase space, 
where complicated self similar fractal structures are present J0|,pp. In fact, 
for such systems, a convex combination of the power law and the exponen- 
tial law decays seems to provide a very good fit of F x {t) and is theoretically 
justified [TTj . 

Our aim is to study the statistics of the first return times for smooth 
perturbations (non necessarily invertible) of the uniform rotations T^q on 
the circle S 1 , 

e^T k , n (9) = 9 + n + ^f(9), (5) 

where 9 G S 1 and k G R + and Q G R + are the parameters of the map. The 
function f{9) is a trigonometric polynomial such that the maps are monotonic 
for k G [0, 1) and non- invertible for k > 1. The details of the dynamics and 
the structure of the bifurcation diagram have been thoroughly studied for 
the sine-circle map given by Tk,n{9) — 9 + Vt + -^sm.2ii9 (see for example 
|14j - [2E])- Other families of the form (4) are used to study the universality 
of the properties found for the sine-circle map. 

The main results of our analysis of the recurrence times for the circle maps 
of the form (5) can be summarized as follows. For k < 1, where the map 
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is invertible and diffeomorphic to a rotation, if the rotation number is dio- 
phantine, three return times are observed. This is in agreement with Slater's 
theorem [T2*] . which we extend from linear rotations to the diffeomorphisms 
of the circle, see section 2. At the critical value k = 1 three return times are 
still observed and the average return times allow an effective reconstruction 
of the invariant measure, see section 3. In the super-critical case (section 4), 
we show the appearance of a continuous spectrum F$(t) for maps that have 
chaotic orbits at k = 1 + e (section 4.1), and for maps T^q for k beyond 
the corresponding accumulation point of period-doubling bifurcations (sec- 
tion 4.3). Also, the properties of the return times could be used to detect 
the existence of attracting periodic orbits (section 4.2). 

Our conclusion is that spectrum reflects, the bifurcations in the topolog- 
ical properties of the dynamics, and is a useful tool to investigate them. 

2 Dynamical properties of the circle map 

Let us briefly recapitulate some of the properties of the circle maps (5). For 
our purposes, we distinguish three regions depending on the parameter k: 
the sub-critical region k < 1, the weakly super-critical region k = 1 + e and 
the strongly super-critical region k » 1 . 

2.1 Sub-critical and critical region 

For k < 1 the circle map T^n is an orientation preserving homeomorphism 
of the circle, and for k < 1 the map is a diffeomorphism. In any case, its 
topological properties are fixed by the rotation number ui defined by 

u = Hm (6) 

where T is the lift of T on the real line, and < ui < 1 for the definition to 
be unique. 
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The map T^u, for k < 1 is conjugate to the linear rotation by the angle 
u(k,Q). The properties of the conjugation = &k,n{9) with the linear 
rotation R w : — > + u depend on the arithmetic properties of u. For 
a generic rotation number $jt,n is a homeomorphism according to Denjoy's 
theorem |14j . 

Furthermore if u is a Diophantine irrational the conjugation $>k,n is a 
diffeomorphism [T^]. In such the topological and metric properties of 

the linear rotation R w extend to the map T^q by using 

Tin = *J$* ° K ° **,n (7) 

Slater's theorem [TS], stating that for any irrational linear rotation, and any 
connected interval there are at most three different return times, one of them 
being the sum of the others, extends from linear rotations to the map T^n- 
Two of the three return times are always the consecutive denominators in 
the continued fraction expansion of the irrational rotation number u, and 
one of the return times is a sum of the other two. Two points 9, T£ n (6) 
in a connected interval A are mapped into G = $fc,fi(0) and of the 

connected interval <&k,n{A) 

A critical map Tk=i,n is only a homemorphism of the circle but is still 
characterized by a unique rotation number. Our calculations indicate that, 
at least, for sufficiently irrational rotation numbers there are still only three 
return times, like in the sub-critical case. 

2.2 Weakly super-critical region 

When k > 1 the map ceases to be invertible and it does not have a unique 
rotation number. For a given (k > 1,0) quasi-periodic, chaotic, at most two 
stable periodic orbits, and orbits asymptotic to the latter can coexist. The 
points in the parameters plane (k, fl) that correspond to sub-critical maps 
with rational P/Q rotation numbers form domains, the P/Q tongues, that 
can be extended above the critical line k = 1 where they start to overlap. 
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The boundaries of the tongues, correspond to tangent bifurcation, and can 
be found by solving for 9 and Q the equations 

T$(0) = + P, T«>(0)>) = 1. (8) 

where T^(6) is the Q-th application of the map, T^(6) its derivative (see 
for example |19j). 

The union of all Pj Q tongues at k — 1 has full measure, and for k = 1 + e 
any map (fc, fi) belongs to intersection of two tongues P\jQ\ and P2/Q2 that 
correspond to Farey neighbours with sufficiently high Qi and Qi- Actually, 
(k, fl) is in the intersection of the tongues that correspond to all rationals 
deeper in the Farey tree and between P\/Qi and PijQ2- To predict the 
dynamics of the map one has to know a very complicated fine structure of 
bifurcations inside the overlapping tongues. The lines in the (k < 1, Q) plane 
that correspond to maps with irrational rotation numbers become, for k > 1, 
domains such that a map in such a domain have at least one orbit with 
the corresponding irrational rotation number, but also has orbits with other 
rotation numbers consistent with the P/Q tongues that overlap, and chaotic 
orbits [201 ■ Weakly super-critical maps of the form (5) with chaotic orbits 
have been used to study the quasi-periodic rout to chaotic dynamics [T5] . 

MM- 

2.3 Strongly super-critical region 

Besides the maps (k, Q) with the small k > 1 that have chaotic orbits, there 
are such Q that maps (k, Q) inside a P/Q tongue will have strongly chaotic 
behaviour only if k > K C (P/Q), where K C (P/Q) is the value of k at which 
the period doubling bifurcations inside P/Q tongue accumulate, which could 
be quite large. For example, for Q = the critical value of k is estimated to 
be k = K c (0) = 4.60366 [15] . We shall see that the properties of return times 
change abruptly at k = K C (P/Q), and can be used to estimate K C (P/Q). 
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3 Return times for sub-critical and critical 
dynamics 



The numerically observed return times for k < 1 and diophantine rotation 
numbers uj are in agreement with a straightforward extension of Slater's theo- 
rem. The same result is found at the critical value k = 1. For generic interval 
and each quadratic irrational uj that we have studied, there are again only 
three return times, one is the sum of the other two and two always coincide 
with the denominators of the corresponding two successive approximants of 
uj, see for example figure la. We are fairly confident that the conclusions 
are valid for any quadratic u and the maps in the class (5), but we can not 
claim anything for homeomorphisms with a generic rotation number (see the 
reference |28j). 

The three return times, and their relative weights depend on the location 
and the size of the interval. However, there is a sequence of intervals, obtained 
by partitioning the circle with the iterates of an initial point 6 that is best 
suited for the analysis of the return times at a point 9q of the map with a 
given irrational rotation number uj. One considers the trajectory formed by 
the qi-i + qi — 1 successive iterates of a point 9 , where Pi/qi are rational 
continued fraction approximates of the uj. 

The qrth and the gj_i-th iterate form the boundary of an interval A^Oq) 
which contains only the initial point 9q and no other points of the considered 
part of the trajectory. The points generated by q i+1 + q { — 1 iterates will 
subdivide the intervals generated by % + q^\ — 1 iterates. One obtains a 
sequence A i+ i C A$ of intervals that converge to the point 9 on the orbit of 
Tk,n- Calculating the return time for such sequence of intervals A4, i — > 00 
is best adopted to the calculation of the return time at the point 9q G D£i A4. 
Suppose that the gj_i-th iterate is to the left of #0 which is to the left of the 
qrth iterate. Then, each of the points 9 e A^Oq) that are on the left side of 
9q have a unique first return time equal to and points 6 e Ai(9 ) to the 
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right of #o also have a unique first return time equal to qi-\. Thus, there are 
only two return times equal to and gj. Furthermore, the union of the 
intervals formed by g^+i + qi — 1 iterates gives a partition V% of the circle, 
and the return times into various intervals of a partition V{ are the same q± 
and qi-i, although the relative weights are obviously different for a nonlinear 
map. However, as 1 — > oo, and for an irrational number uj = [a, a, a, . . .] with 
a constant tail of the continued fraction expansion, the relative weights of 
the two return times become independent of i, which implies the existence 
of a point spectrum independent of 6 (see Appendix) 

F(t)=w-8(t-t-)+w+5(t-t+), (9) 

where 

t- = lim — — = — - — t—, t+ = h m — ^— — = wt-, (10) 

are the renormalized return times in the limit and are denominators of the 
continued fraction approximants of uj. Since the map is smoothly conjugated 
for k < 1 and diophantine u to a linear rotation, the above properties follow 
from a proposition we prove in this case in the Appendix. 

Approximation of the measure: The distribution of the return times in 
various intervals can be illustrated using the mean return time. In fact, for 
any homeomorphism of the circle there is a unique invariant ergodic mea- 
sure which is absolutely continuous with respect to the Lebesgue measure 
on the circle, and the density of this measure can be obtained using return 
times, as is indicated by the Kac lemma. In figures lb,c,d we plot density 
of a coarse-grained mean return time, i.e. the ratio between uniform aver- 
age of the return times into an interval of a partition of the circle divided 
by the Lebesque measure of the interval, for a sub-critical and the critical 
sine-circle map. For a sufficiently fine partition this quantity illustrates the 
density of the unique invariant ergodic measure for the considered map. The 
main computational cost of this method is due to the computational time, 
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complementary to the perturbative method where the main requirement is 
for a sufficient storage space. In order to achieve an resolution of 10 _m by the 
return times calculations one needs roughly 10 m x K10 m x 10 2m iterations of 
the map. Here, 10 m is the number of intervals of the partition, K10 m with 
K < 10 is roughly the average return time and 10 2m is the number of points 
in each of the intervals. The same space resolution is archived by using 10 m 
Fourier components of the conjugation function in the perturbation method. 
For 1 percent accuracy, i.e. m = 2, both methods can be easily implemented, 
but 10~ 3 accuracy represents a more challenging task for both the methods. 

For k < 1 the density of the invariant measure is a smooth function (see 
figure lb,c), and for the critical maps the density becomes singular (figure 
Id) |2I3>[ZS]- Fractal properties of the ergodic measures for the critical circle 
maps (5) have been studied. There is strong evidence that the class 
of critical maps with the same fractal spectrum of the invariant measure are 
characterized only by the rotation number (actually probably only by the tail 
in it's continued fraction expansion and by the type of the singularity 
that induces the critical behaviour. 

4 Return times in the super-critical cases 

A super-critical map could have chaotic orbits, at most two stable attracting 
orbits and orbits asymptotic to these. Our numerical computations sup- 
port a conclusion that for any Q there is sufficiently large k, such that the 
distribution of the first return times is given by exponentially fast decay. Fur- 
thermore, there are maps T^q which have chaotic orbits for k — 1 + e with 
arbitrary small but non-zero e, and maps Tk,n that show chaotic behaviour 
only for sufficiently large k > K C (P/Q). This is manifested by two different 
roots to the exponential decay of the distribution of the first return times. 

The following three typical situations can be clearly distinguished by 
studying the properties of the return times. 
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4.1 Quasi-periodic route 

Consider, first, a weakly-super-critical map (k, Q) with k — 1 + e, where e is 
arbitrary small but non-zero, and with Q such that there is an orbit with an 
irrational rotation number. For example, suppose that the rotation number 
uj(0q) of the orbit through point O = is numerically equal to the golden 
mean 7. We can use the return times into a sequence of nested intervals that 
shrinks to 9 , or analogous sequences that shrink on other points to study 
the dynamics. 

Results of such analysis are shown in figures 2a,b,c,d which illustrate the 
dynamics of the same map (ko,Qo) but as revealed by the statistics of the 
return times into a sequence of nested intervals of decreasing size. Here 
ko = 1.01 and Q = 0.606494989, leading to an orbit with uo whose continued 
fraction expansion is u = [0, 1, 1,1,1, 1, 1,1,1, 1, 1, 1, 2, 1, . . .], and which has 
denominators of the convergents = 2, 3, 5, 8, 13, 21, 34, 55, 89, 144, 233 . . .. 

The fact that the map is not topologically equivalent to a uniform rotation 
is manifested already at the resolution given by the interval of finite size. 
For example, for the interval (0, 0.0015) (fig. 2a) there are four return times 
55, 89, 144, 199, where the fourth 199 is not a denominator of any of the 
convergents to u. However, it is the denominator of the Farey neighbor 
123/199 of the approximant 89/144. This signals that, for k = k the 89/144- 
tongue and the 123/199-tongue have common part of their interiors, and that 
the point (ko, Qq) belongs to both tongues. The interval that will detect the 
existence of two intersecting tongues at (ko, Qq) must be smaller than the 
distance between 9 and its 144-th iterate. The first return times into a 
larger interval are at most 144 for all its points, i.e. with the resolution 
weaker than the points of the first 144 iterates of 9 the map looks as a 
smooth rotation. 

Further analysis of the return times for the same map (ko, ^o) but into 
smaller intervals reveals intersections of tongues at the deeper levels of the 
Farey tree between 89/144 and 123/199. This is illustrated in figures 2b, c 
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which show the return times 89, 144, 199, 233, 288, 343 into the interval (0, 0.001), 
and a very large number of return times into the interval (0,0.0001). The 
statistics of the return times into A = (0, 0.0001) is illustrated in figure 2d, by 
plotting the logarithm of the probability density In F do (t) of the return time 
larger that t — rj < ta > versus the normalized time t. The distribution 
is given by exponential decay with an exponent that is numerically close to 
I. Thus, on the sufficiently small scale, the map has the distribution of the 
return times characteristic of strongly chaotic systems. 

Increasing k, and changing Q so that there is always an orbit with rota- 
tion number 7, moves the point (k, Q) into the domain were the intersections 
of tongues, and the chaotic behaviour, can be detected using larger intervals. 
Figures 3a,b,c, show the effects of more intersected tongues on the distribu- 
tion of return times into the same interval A = (0, 0.0015). For k = 1.2, the 
return times and their relative weights are such that the distribution Fg (t) 
is given by exponential decay with the exponent that is again numerically 
equal to 1. In fact, for any k > 1 one can find pairs (k, Q) which imply non- 
periodic orbit through the point 9 , and the the density F 9o (t) for such map 
(k, f2) is F do (t) = exp(— t). This is illustrated in figure 4a for k = 1.3; 3; 4; 5 
and k = 6.14 and the corresponding f2. For all these maps and for all tested 
9 the distribution F e (t) is always Fg (t) = exp(— t). 

Numerical evidence, presented in figures 2 and 3, suggests the following 
conjecture: Suppose that a point 9q lies on a quasi-periodic or on a chaotic 
orbit of a map {k = 1 + e, Q) , for arbitrary e 7^ and consider a sequence Ai 
of nested intervals containing 9 whose length approaches zero by increasing 
i. The number of different return times also increases with i, and asymptoti- 
cally, as Ai shrink to 9, the distribution of return times becomes continuous. 
Furthermore, the density of probability F e (t) with respect to the Lebesgue 
measure on the circle of the normalized first return time larger than t is given 
by exponential decay. 



11 



4.2 (Bi)-stability 

Possible existence of attracting periodic orbits can be detected by studying 
the return times into various intervals of a single partition of the circle. 
Although there could be no bounded invariant density in this case the return 
times are still well defined. The return times into different intervals depend 
on whether the stable orbit have points in the considered interval or not. 
If there is an attracting periodic orbit with no points in the interval than, 
obviously, there are points in the interval that will never come back, i.e. with 
the first return time r = oo, indicating the existence of the attracting orbit. 
In this case there will be only those return times that correspond to orbits 
that re- visit the interval at least once before being attracted to the attracting 
periodic orbit. On the other hand, for examples of chaotic maps, illustrated 
in the figure 4a, the return time statistics for increasingly fine partitions 
confirms that these maps have no attracting periodic orbits. 

4.3 Period-doubling route 

In order to study how the period-doubling route to an ultimately chaotic map 
(k,Q = P/Q) is manifested in the properties of the return times, consider 
the maps (k, Q = 0) for various k. Results are illustrated by various curves in 
figure 4. The return times into intervals at 8q = for any k G (1, K C (Q = 0) 
show the existence of stable periodic orbits, as described in the previous sub- 
section. Suddenly, at the accumulation points of the period-doubling cascade 
k = K C (Q = 0), the distribution of the return times becomes continuous. For 
such critical k, the distribution Fg (t) at the point 8q = is given by the ex- 
ponential decay aexp(— at), with a < 1 and a 0.8 ^ 1. Furthermore, 
still for k = K C (Q = 0) the return times into intervals at 6q of finite size 
are given by double exponential curves, which converge to the single expo- 
nential exp(— at) as the intervals shrink to 9 . For example, for the interval 
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A = (0,0.01) the distribution F A (t) is well approximated by 

F A (t) = 0.15exp(-0.257t) + 0.7exp(-1.8t) (11) 

represented by the dash-ed line in figure 4b. 

Other curves in figure 4 represent the distributions of the return times 
for examples of strongly super-critical maps, and for intervals of decreasing 
size around different points. The curves with the unique slope in figure 4b 
are the distributions for maps in the tongues with Q = for k beyond the 
accumulation of period doublings, i.e. for k > 4.604 and with Q = 1/2 for 
k > 1.978. In all the cases, the distribution is given by exponential decay 
that converges to F 6o (t) = aexp(— at), where a ^ 1 and 9 is on an unstable 
periodic orbit. Non-linear curves in figure 4b represent the distributions 
for the map at the accumulation of period doublings and on the indicated 
intervals of decreasing size. The curves in 4a represent F e (t) for examples 
of strongly super-critical maps at 9 not on a periodic point, when F do (t) = 
exp(— at) with a m 1.0. Shown are examples with k = 1.3,3,4,5,6.14 and 
the corresponding Q as explained earlier, and also examples of maps with 
il = 0; 1/2 and large k. All this is consistent with the statistics of the return 
times for other examples of strongly chaotic maps. 

5 Summary and conclusions 

We have analyzed the circle map numerically for a wide range of values in 
the parameter space {k, Q) and different sets of initial conditions. The results 
can be summarized as follows. 

Sub-critical and critical region: For k < 1 three return times are observed 
and this is theoretically and numerically justified. For diophantine rotation 
numbers oo this result can be presented as a corollary of Slater's theorem, 
since the map is diffeomorphic to a linear rotation, to which such a theorem 
applies. For the case of a special sequence of nested intervals including a 
given point we provide a very simple proof of Slater's theorem, showing 
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that there are only two return times and proving the existence of a limit 
point spectrum F(t). The critical dynamics k — 1 is also clear since three 
return times for a generic interval and for each quadratic irrational rotation 
number are observed. For k < 1 and u diophantine a piece-wise constant 
approximation to the invariant measure is obtained from the average return 
times from a uniform partition of the circle with a very simple procedure and 
an accuracy comparable to other methods. 

Super-critical region: The dynamics of a map (k, fl) in the weak super- 
critical case k = 1 + e is dictated by the P/Q tongues with a nonempty 
intersection that contains the point (k,fl). This is manifested, and could be 
detected, in the distribution of the first return times by appearance of more 
than three return times, which correspond to the rationals on the Farey 
tree in-between the the two major overlapping tongues that contain (k,fl). 
In the case that the interval contains a point on a non-periodic orbit than 
there is a sub-interval such that the distribution, with respect to the uniform 
distribution of initial points, of the return times into this sub-interval is 
typical for strongly chaotic systems, i.e. the exponential decay with exponent 
equal to unity. 

In the intermediate region two return times or the exponential-one spectrum 
are observed depending on the existence of attracting periodic orbits. The 
way -Fa(0 — ► e~* when the size of A containing 9 approaches depends on 
(k, fl, 6q) and it is convenient to distinguish various routes. 

Quasi -periodic route: For any non-periodic point 9o of the map (k, fi) 
with k > 1 + e where e > is arbitrary small, the spectrum is exponential- 
one. Choosing a finite interval A the spectrum F^i) appears as continuous 
for a sufficiently small A. 

Bi-stability: For k = 1 + e there are values of (k, fl) such that the map has 
attracting periodic orbits. For any interval A not intersecting one of these 
orbits one of the return times is r = oo, since many points do not return 
being attracted by the periodic orbit. For 1 < k < 2 the attractive periodic 



14 



orbits are present for most values of f2. 

Period-doubling route: For any rational Q = p/q there is a critical 
k c = k c (p/q) corresponding to the accumulation of period-doublings in every 
tongue. The transition from (bi)-stability to chaoticity is manifested abruptly 
in the spectrum F^t). For k = k c the spectrum F^if) is continuous and there 
are intervals such that it can be fitted with a double exponential: in the limit 
f^(A) — > the spectrum becomes exponential. 

Strongly super-critical region For k 3> 1 the map is chaotic and the 
periodic orbits are unstable. The spectrum F(t) is exponential-one except 
for a set of points of measure zero corresponding to unstable periodic or- 
bits where the exponential decay speed is different from 1. The results for 
the super-critical dynamics indicate that the analysis of the return times 
spectra in the super-critical case could be a useful tool for a better under- 
standing of the transition from the weakly to the strongly chaotic regime. 
It would be interesting to follow in details the pattern of intersections of 
tongues and period doubling bifurcations inside each tongue, leading to the 
strongly super-critical case, by using the return times spectra and the way 
they are approached when a sequence of nested intervals squeezing to a point 
is considered. To conclude the computation of the return times spectrum is a 
simple and effective way to explore a dynamical system and its bifurcations. 
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6 Appendix 

Letting be the linear map conjugated to and Oo = $ -1 (#o) be the 
image of the initial point do we consider a partition of the circle -R"@o to 
which corresponds another partition T£ n (6 ). The order in these partitions 
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is preserved since the maps are diffeomorphic for k < 1 and u diophantine 
Let the continued fraction expansion of uj be given by 

u=[a 1 ,a 2 ,...] = — (12) 

ai H 



a 2 + ' ■ ■ 

and Pi/qi = [ai, . . . ,cii] be the corresponding rational approximations of 
order %. The odd and even approximants are upper and lower bounds to uj, 
converging monotonically to it. The following recurrence relations hold 

Pi = a,iPi-i + Pi-2 qi = aiqi-i + g;_ 2 (13) 

and the following inequalities hold 

1 11 1 . . 

< p^ - uiqi-x < — — ■ <qiUJ-pi< (14) 



ft-i + qi " Qi Qi + Qi+i ' ' ' Qi+i 

Denoting the linear map iterates of a point O by 

eg* = R q je = e + g^mod 1 = e + qi u-p l (15) 

the odd and even sequences 6o 2l_1 and Q q 2i converge monotonically from 
below and from above to @ . The intervals A(©o) = [@o ~\ @o ] if * is even 
and Ai(Q ) = [0g% 0q _1 ] if % is odd, form a nested sequence of intervals 
Ai D . . . Ai D A i+1 . . . squeezing to O as i — > oo and from the previous 
inequalities the following inclusions hold. 



©o ©0 + ! 



+ Qi ' Qi + Qi+i 



c M9 ) c 



©o - -, ©o + ' 



Qi qi+i 



(16) 



The intervals B { = ^ k } n (Ai) = [T l k J(9 ),Tl n (9 )} enjoy the same proper- 
ties for any diophantine u since the map is orientation preserving and &k,n 
is a diffeomorphism. The sequence Bi is a nested monotonic sequence of 
intervals approaching #o exponentially fast. According to Kac's lemma the 
average return times for the intervals Bi and Ai are given by the inverse of 
their length, which increases to exponentially fast with i. 
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According to Slater's theorem for a generic interval and a linear map with 
irrational u>, there are three return times, the last one being the sum of the 
first two. The sequences of nested intervals [©o~\ @o] and [Go, ©o] enjoy 
this property. For the intervals Ai defined above the return times are only 
two and we give a sketch of the proof since it quite simple. 

Proposition: The return times in the interval Aj(0 o ) for the linear map 
with an irrational u are two. If i is even (odd) and for Q G A^Qq) the return 
time is if < 0o and qi-i (<&) if > O . 

Proof The proofs for even or odd % are analogous so we consider only the 
even i. In order to prove the results we suppose O is sufficiently far from 
the identified ends and 1 that the order preserving relation becomes the 
inequality between real numbers. For initial conditions near 0, the torus 
defined as the interval [—1/2, 1/2] with identified ends should be considered 
to have the same correspondence. The Lagrange's theorem states that in the 
interval 1 < n < qi + q^i — 1 the minimum of \nuj — m\, where m = [nuj] is 
the integer part, is reached for n = qt. As a consequence, the two points 0q 
closest to O for n < + g^-i — 1 correspond to n — <&_i, qi- We consider the 
sequence of nested intervals Ai(Q ) = [0o ~\ 0o9i] an d for a fixed i a point 
G Ai(Q ) we denote by n = i?"0 the points of its orbits and examine 
two possible cases. 

Case 1 If < 0o then: 0" £ At for n < q u 9i G A, t . 

To prove the first point we notice that 

e* - e qi = O - > o (17) 

and 

QQi _ Q«-l = Q 9i _ Q«-i _ ( «-l _ 0ffi-i) = ^ A .^ _ ( @0 _ @) > (18) 

since both O , belong to A4. To prove the second property we notice 
that the points 0„ closest to for n < qi are met for n = qi-2, q%-i and 
Qqi-i < < Qn-2 anc l we s how that both are out of A^. 
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eg*- 1 - e^- 1 = e - e > o (19) 

To prove Q q ^ 2 > Q q l we show that Q q ^ 2 - > ^(A)- Indeed /i L (A) = 
ujqi -Pi- ujq-i-i - Pi-i = uj(qi - - (pi - ft_i) 

= uqi_ 2 ~ Pi-2 + {fli ~ 1) (w qi-i - Pi-i) (20) 

Consequently we obtain 

® q '- 2 - © = wqt-2 - = A*l(A) + (at - 1) (Pi_i - u > /x L (A) (21) 

the equal sign holding for the golden mean. 

Case 2 If > O then: 9™ £ A { for n < G A*- 

The fist point is proved observing that 

e9l -i _ e «-i = @ _ @ o > o (22) 

and 

09, _ 0<7,-i = Q <n _ Qgi + 9l _ 9i _i = ^(^ _ (0 _ @0 ) > o (23) 

Concerning the last point we notice that we already proved (two equations 
above) that Q qi -' 2 — > /i L (Aj) for any e i,. We prove also that — 
Q q '- S > fi{Ai). Indeed from - Q q '-* - /x(A;) = ft__ 3 - c^_ 3 - (wft - ft) + 
(ujqi-i - Pi-i) = 

= Pi- a,i_i p^ 2 -uj(qi- ai_igj_ 2 ) > Pi ~ Pi-2 ~ w( % - ft-2) 

= aj(ft-2 _ ^<7i-i) > 0, where we have used twice the 
recurrence relations for the pi and taking into account that gj_2 — pi-2 > 
0, and that > 1 since we have assumed u to be irrational. 

Existence of the spectrum: The return times in the intervals Aj are 
given by q { if G A~ = A0 G [©o" 1 ,^] and if G A+ = [0 o ,©o]- 
As a consequence in we consider the return times in Ai the relative measures 
of points that come back to Aj — and to Af will be given by the Lebesgue 
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measure of these intervals normalized to the Lebesgue measure of A^. As a 
consequence the weights W-, w + for the return times qi and q^\ are given 

by 

fi L (Ai) l + r + fi L (A) l + r- 1 /ii(4 r ) 

We notice that these weights for a quadratic irrational are independent of %. 
Indeed in this case pi = g £ _x and letting u = [a, a, a, a, . . .] 

r = "ft-*- 1 = -J— (25) 
g,_ 2 - to a + 

Indeed taking into account that cj(a + cj) = 1, the above relation is verified 
if 9j — (a + a>)gi_i — — 2 + cugi — 1 = ^ — agj_i — gj_ 2 , which is the case 
due to the recurrence for q { . As a consequence the average time is 

/ \ , (a + u) qj + 

(taJ = + w+gj-i = — ■ (26) 

1 + a + uj 

The normalized times in the limit % — > oo become 
i_ = hm = (w + + ujw + ) 



I \ V — )— I ^aj _j_ y 

(taJ ' a + 2cu 



t+ = hm - — - = u>_ + a;u» + = —cut- 27 

i^oo ( TA .) v +J a + 2u v 7 

As a consequence the limit point spectrum exists and is given by 

F(t) = w-5(t-t_)+w+5(t-t+) (28) 

We believe that the spectrum exists for any irrational uj. The existence of a 
spectrum for a generic nested sequence of intervals, squeezing to 9 , remains 
an open question. 
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FIGURE CAPTIONS 

Figure la,b,c,d:Illustrate the sub-critical and critical cases with u = 7 
:(a)Three return times for k — 1; ft = 0.606661; A = (0, 0.0015), and densities 
of the unique invariant ergodic measure for (b) k = 0.75, fl = 0.61088669 (c) 
k = 0.9, D, = 0.6083938 and (d)k = l,Q = 0.606661. 

Figure 2a,b,c,d: Illustrate the distribution of the return times into a se- 
quence of sub- intervals for the fixed weakly super-critical map k = 1.01; Q = 
0.606494989. 

Figures 3a,b,c,d:Illustrate the return times into the interval (0, 0.0015) 
for weakly super-critical maps with (k, ft) = (a) (1.015,0.6063931) (b) (1.1,0.6054765) 
(c) (1.2,0.6047099) (d) shows kiF e (t).vs.t for k = 1.2, Q = 0.6047099, with 
the slope 0.9996. 

Figures 4a,b: Strongly super-critical dynamics: (a) ]n.Fe(t).ys.t with 
slopes 1.0 for (k,Q) = (1.3, 0.606187), (3, 0.51738887), (4, 0.36176849), (5, 0.20197433), 
(6.14, 0.7184975) and for (k, Q) = (4.605, 0), (211, 0),(1.98, 1/2), (2, 1/2); 

b)lnF A (t).vsi with slope ^ 1, for (k, ft) = (4.60366,0) when A = (0.01); 
(0,0.0005) and A = (0,0.0001), with approximation (14) (dash-ed line). 
Other curves represent lnF^^(t).vs.t at unstable periodic points for other 
examples with ft — 0; 1/2. The maximal slope is 0.85. 
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